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The free energies of static quarks and the Debye screening masses in the quark gluon 
plasma are studied using Polyakov-line correlation functions in lattice QCD adopting the 
i i. fixed-scale approach in which temperature is varied without changing the spatial volume and 

the renormalization factors. We calculate static-quark free energies in various color channels 
in the high temperature phase up to about 3.5 times the (pseudo-) critical temperature, per- 
forming lattice simulations both in quenched and 2 + 1 flavor QCD. For the quenched simula- 

■ tions, we adopt the plaquette gauge action on anisotropic 20 3 x Nt lattices with Nt = 8-26 at 
the renormalized anisotropy a s /a t — 4. For 2 + 1 flavor QCD, we adopt the renormalization- 
group improved Iwasaki gluon action and the non-perturbatively 0(a)-improved Wilson 
quark action on isotropic 32 3 x Nt lattices with Nt = 4-12 at mps/mv = 0.63 (0.74) for 

■ the light (strange) flavors. We find that the color-singlet free energies at high temperatures 
' converge to the zero-temperature static-quark potential evaluated from the Wilson-loop at 
_ short distances. This is in accordance with the theoretical expectation that the short dis- 

■ tance physics is insensitive to the temperature. At long distances, the free energies approach 
• • ■ to twice the single-quark free energies, implying that the interaction between static quarks 

is fully screened. We find that the static-quark free energies for various color channels turn 
k>( , out to be well described by the screened Coulomb form, and the color-channel dependence 

of the inter-quark interaction can be described by the kinetic Casimir factor inspired from 
the lowest order perturbation theory. We also discuss comparison with a prediction of the 
thermal perturbation theory and flavor dependence of the screening masses. 

§1. Introduction 

At temperature around 150-200 MeV, the hadronic matter transforms into the 
quark-gluon plasma (QGP) by the QCD phase transition. Recent experiments of rel- 
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ativistic heavy-ion collisions at RHIC and LHC make it possible to directly explore 
the properties of QGP. 1 ' Because the interaction is strong, lattice QCD simulations 
are indispensable in theoretical investigations of the QGP. So far, most lattice sim- 
ulations at finite temperature have been performed by using staggered-type lattice 
quark actions^ iGDSli©© with which the computational cost is less demanding. How- 
ever, the fourth-root procedure which is required to make a realistic quarks by the 
staggered-type actions is theoretically not fully established yet. In this situation, it 
is of great importance to check the results by other actions such as the Wilson-type 
quark actions, which have solid theoretical basis. 

A series of systematic studies of finite-temperature QCD with the 0(a)-improved 
Wilson quark action was made by the CP-PACS Collaboration^® With two flavors 
of degenerate dynamical quarks, the phase structure, the pseudo-critical tempera- 
ture, the 0(4) critical scaling around the chiral transition, as well as the equation 
of state have been investigated. The WHOT-QCD Collaboration is extending the 
project in various aspects — to more detailed properties of QGP such as the free 
energies of static quarks and the electric and magnetic screening masses GDflEJiinj to 
2 + 1 flavor QCD including the dynamical strange quark^J'^ and to finite densi- 
ties.^'^ In this paper, we study static-quark free energies in various color channels 
as well as their Debye screening masses in QCD with 2 + 1 flavors of improved Wilson 
quarks. 

The free energies for static quark-antiquark and two static quarks characterize 
inter-quark interactions in QGP, and their Debye screening masses describe the ther- 
mal fluctuation of quarks and gluons in QGP. In a phenomenological model, they 
play important roles in the fate of heavy-quark bound states such as J/ip and T in 
QGP created at the center of heavy-ion collisions at RHIC and LHCP^ On the 
lattice, the static-quark free energies are extracted from Polyakov-line correlation 
functions. Studies in quenched QQjj^^^^ and in full QCD with staggered- 
type quark actions^ or with Wilson-type quark actions^ MXM h ave been 
made, and comparisons with analytic studies^ '^ have been attempted. The stud- 
ies with Wilson-type quarks were limited to the case of two-flavor QCD. We now 
extend the study to 2 + 1 flavors. 

To reduce the high computational cost needed to simulate Wilson- type quarks, 
we adopt the fixed-scale approach J^ 1 ED [ n which the temperature T = 1/aNt 
is varied by changing temporal lattice size Nt with fixed values of the coupling 
parameters (and thus with a fixed value of the lattice spacing a). Because the 
renormalization depends only on the coupling parameters, we can use a common 
zero-temperature simulation to obtain non-perturbative renormalization factors to 
renormalize observables at all temperatures. We may even borrow existing config- 
urations generated for zero-temperature studies. We can thus reduce the costs for 
the zero-temperature simulations, which is quite demanding in finite-temperature 
studies. 

In the fixed-scale approach, we vary T without changing the system volume and 
lattice cut-off effects. We can thus study the pure effects of T. We show later that 
this property of the fixed-scale approach benefits us also in a study of static-quark 
free energies. On the other hand, results of the fixed-scale approach at very high 
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temperatures, typically at T ^ l/3a, suffer from lattice artifacts due to coarseness 
of the lattice to resolve thermal fluctuations (see discussion in Sec. 14. ip . We discuss 
this point later in this paper. 

We first make test calculations of the static-quark free energies with the fixed- 
scale approach in quenched QCD. We then study the case of 2 + 1 flavor QCD 
adopting a non-perturbatively improved Wilson quark action and a renormalization- 
group improved Iwasaki gluon action at mps/my = 0.63 (0.74) for the light (strange) 
flavor, where mps (my) is pseudo-scalar (vector) meson mass at T = 0. Comparison 
to the quenched results is made to identify the effects of dynamical quarks. The 
lattice actions and the simulation parameters are summarized in Sec. [2J Results of 
the free energies of static quarks in various color channels are given in Sec. We 
show that the color-singlet static-quark free energies at high temperatures universally 
converge to the zero-temperature static-quark potential evaluated from the Wilson- 
loop at short distances, whereas, at long distances, they approach to twice the single- 
quark free energy. These are in accordance with theoretical expectations from the 
asymptotic freedom and the screening of colors in the high temperature deconfined 
phase. Thanks to the fixed-scale approach, we directly confirm these properties. The 
color-channel dependence of the free energies is also discussed. In Sec.[U the effective 
coupling a c fj(T) and the Debye screening mass TO£>(T) are calculated by fitting the 
free energies with a screened Coulomb form. We discuss the effects of dynamical 
quarks and compare the results with a prediction of the thermal perturbation theory. 
A summary is given in Sec. In Appendix [Aj we discuss systematic errors in the 
Debye screening masses. 



§2. Simulation parameters 



2.1. Quenched QCD 

In quenched QCD, we adopt configurations generated in Ref. 29) as the series 
"a2" on anisotropic 20 3 x Nt lattices at the renormalized anisotropy a s /at — 4, 
where a s and at are the lattice spacings in the spatial and temporal directions, 
respectively. The anisotropic lattice enables us to study temperature dependence 
with fine resolution in the fixed-scale approach since the resolution is proportional 
to l/(Ntat). This will be useful to compare the screening mass in the fixed-scale 
approach with those in the previous fixed- Nt approach in Sec. 14.11 In Ref. |29|) . 
results obtained on isotropic lattice with a similar spatial lattice spacing are shown 
to be well consistent with those obtained on the anisotropic lattice. The gauge action 
is the standard plaquette action with anisotropy, 



^ReTr[/i 4 (x) 
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where U fJiU (x) is the plaquette in the [iv plane and (3 and £o are the bare lattice gauge 
coupling and bare anisotropy parameters. Tuning to the renormalized anisotropy 
£ = 4 as well as zero-temperature simulation on an anisotropic 20 3 x 160 lattice was 
done in Ref. GSty . 
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At our simulation point ft = 6.1, the spatial cutoff scale is a s ~ 0.097 fm and 
the spatial length of the lattice is about 1.9 fm, where the lattice scale is set by 
the Sommer scale tq = 0.5 fmP" We use the configurations at Nt = 8, 10, • • • , 26 
in the high temperature phase. This range of Nt corresponds to the temperature 
range 1.08 & T/T c £ 3.50, as summarized in Table HI where T c 290 MeV is the 
critical temperature of the pure gauge theory. Up to a few million configurations 
were generated at each Nt using the pseudo-heat-bath algorithm. Although some 
configurations in the low temperature phase are also available, it turned out that the 
statistics is not high enough to extract static-quark free energies. We thus confine 
ourselves in the high temperature phase. Results for the equation of state with this 
action are given in Ref. | 



2.2. 2 + 1 flavor QCD with Wilson quarks 

For the simulations in 2 + 1 flavor QCD, we employ the renormalization-group 
improved Iwasaki gluon actiorPD coupled with 2 + 1 flavors of the 0(a)-improved 
Wilson quarks! 3 -^* 

S = Sg + Sq, 

(4 4 

f=u,d,s x,y 

where c\ = —0.331, Co = 1 — 8ci, and W™^™ is the n x m rectangular Wilson loop. 
The quark kernel is given by 

M l,y = S x,y ~ Kf {( X ~ lv) U x,»8 X +(i,y + (1 + TmJ^m^W+a} 
— d x ,y Kf CSW £ Viu>Fx,iu>, 

where Kf is the hopping parameter for flavor / and F x ^ v is the lattice field strength 
given by the standard clover-shaped combination of gauge links. The clover coeffi- 
cient csw is determined non-perturbatively as a function of 

On isotropic 32 3 x Nt lattices with Nt = 4, 6, 8, 10, and 12, we perform simu- 
lations at ft = 2.05, = (0.1356,0.1351) and csw = 1-628, which correspond 
to mps/my = 0.6337(38) for light quarks and mpg/my = 0.7377(28) for the strange 
quark. 341 The lattice spacing is a ~ 0.071 fm from the static-quark potential at 
T = (see Sec. EJ). The range N t = 4-12 correspond to T ~ 230-700 MeVpJ The 
results of Polyakov loop and its susceptibility as well as the equation of state suggest 
that the pseudo-critical temperature is around 190 MeV at our simulation pointP^ 

After thermalization, we generate configurations by the HMC+RHMC algo- 
rithm. The Metropolis test is performed every 0.5 trajectories for finite-temperature 
simulations. We measure observables at every 5 trajectories. The number of tra- 
jectories and corresponding temperatures for each Nt are summarized in Table HU 
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Table I. Temperature with respect to T c ~ 290 MeV for each Nt in the quenched QCD simulations, 
where T c is the critical temperature of the pure gauge theory. For each temperature, we use 
100-400 configurations to measure observables. L is the Polyakov loop. 



N t T/T c (L) 



26 


1 


08 


0.0542(5) 


24 


1 


17 


0.0784(5) 


22 


1 


27 


0.1051(3) 


20 


1 


40 


0.1371(4) 


18 


1 


56 


0.1767(4) 


16 


1 


75 


0.2249(4) 


14 


2 


00 


0.2847(3) 


12 


2 


33 


0.3577(5) 


10 


2 


80 


0.4462(2) 


8 


3 


50 


0.5507(3) 



Table II. Number of trajectories and temperatures for each Nt in 2 + 1 flavor QCD simulations. 
Temperatures are estimated using the Sommer scale ro = 0.5 fm. L is the Polyakov loop. 



N t Traj. T (MeV) (L) 



12 


6460 


232 


0.00491( 6) 


10 


3935 


279 


0.01470(11) 


8 


2565 


348 


0.04072(12) 


6 


2495 


465 


0.10981(11) 


4 


3510 


697 


0.29185( 7) 



The statistical errors are estimated by the jackknife method with the bin-size of 20 
trajectories.^ 

§3. Free energy of static quarks 

We now study the free energy of static quarks at zero and finite temperatures. 

3.1. Static-quark potential at T = 

At T = 0, interaction between static quark (Q) and antiquark (Q) can be de- 
scribed by the static-quark potential evaluated from the Wilson loop operator; 

V(r) = - lim [r 1 In (wL xi )] , (3-1) 
I— >oo L \ /J 

where r is the distance between Q and Q. It is expected that the static-quark poten- 
tial shows the Coulomb-type behavior at short distances due to perturbative gluon 
exchange, while it takes the linear behavior at long distances due to the confinement. 
Therefore, its phenomenological form may be written as, 

V(r) = -— + ar + V . (3-2) 

r 
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Fig. 1. Free energies of static quarks in color-singlet channel of quenched QCD at finite temperature. 
The dashed gray curve shows the static-quark potential V(r) at T — calculated in Ref. l30p . The 
arrows on the right side denote twice the single-quark free energy defined by 2Fq — —2T ln{Tri7). 
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Fig. 2. The same as Fig. [T] but for 2 + f flavor QCD. The static-quark potential V(r) at T = 
has been calculated by the CP-PACS and JLQCD Collaborations.^ The fit results of V(r) by 
the Coulomb + linear form are shown by the dashed gray curve. The arrows on the right side 
denote 2F Q = -2Tln(TrJ?}. 



The Sommer scale r$ is defined as r 2 dV/dr\ r=To = 1.65. We set the scale by a 
phenomenological value ro = 0.5 fm. 

In quenched QCD, a fit to the simulation results obtained on the anisotropic 
20 3 x 160 lattice with Eq. fl32J) givesP 

(a a s /a t , aa s a t , V a t , r /a s ) = (0.068(2), 0.0132(1), 0.170(1), 5.140(32)), 

where the fitting range was chosen to be 1.73 < r/a s < 10. Setting ro = 0.5 fm gives 
l/a s = 2.030(13) GeV (a s ~ 0.097 fm) and ~ 0.47 GeV. The fit result of the 
static-quark potential is shown by the dashed gray curve in Fig. [TJ 

For 2 + 1 flavor QCD, the static-quark potential was measured by the CP- 
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PACS and JLQCD Collaborations with the same simulation parameters on the lattice 
N* xN t = 28 3 x 56.® Results are shown in Fig. [2] by the "x" symbol. Because the 
string-breaking effect is not manifest in the range of r we study, we fit the data in 
the fit range of 2-4 < r/a < 14 with Eq. (|3 2[) to obtain 

(a , era 2 , V a, r /a) = (0.44(3), 0.024(1), 0.801(5), 7.06(3)). 

We find l/o = 2.79(1) GeV (a ~ 0.071 fm) and y/a ~ 0.43 GeV. The fit result is 
shown in Fig. [2] by the dashed gray curve. We find that the phenomenological form 
reproduces the lattice data well. 

3.2. Static-quark free energies in the high-temperature phase 
At T > 0, we calculate the Polyakov loop defined as 

N t 

L = Trfi(x) = Tr JJt/^, (3-3) 

T=l 

where i/( X T ) 4 is the link variable in the temporal direction. Results of the Polyakov 
loop expectation value in quenched*^! and 2 + 1 flavor QCD are summarized in Table U 
andHU respectively. The Polyakov loop is related to a single-quark free energy defined 
as Fq = — Tln(L), which expresses the free energy of a static quark in the thermal 
medium. 

The inter-quark interaction at T > may be characterized by the free energy of 
a static quark and an antiquark separated by a distance r, i.e. by the Polyakov-line 
correlation function. The correlation function between Q and Q can be decomposed 
into irreducible representations of color singlet and octet in the color space, i.e. 
3 ® 3* = 1 © 8, whereas that between Q and Q can be decomposed into color 
antitriplet and octet, i.e. 3 © 3 = 3* © 6. Under the Coulomb gauge fixing, these 
correlation functions are given a s,®'^ 



e -^(r,T)/T = I (Tl . 

3 



fl(x)flt( y ) \ (3-4) 



e" F8(r ' T)/T = I (Tr[fl(x)] Tr[flt(y)] ) - 1 (tt [rt(x)flt( y )] ) , (3-5) 
e -f<r,T)/r = J_ <Tr[fl(x)] Tr[fl(y)]> + ^ (Tr[fi(x)fi(y)]> , (3-6) 
e -F3>,T)/T = 1 (Tr[ ^ (x)] Tr ^ (y)]) _ 1 (Tr[^(xMy)]> , (3-7) 

where r = |x — y|. 

Figures Q] and [5] shows the singlet free energy F 1 (r, T) at different tempera- 
tures in quenched and in 2 + 1 flavor QCD, respectively. We find that, i ?1 (r, T) 
at all temperatures converges to the zero-temperature static-quark potential V(r) 
at short distances. This is in accordance with the theoretical expectation that the 



*' In the quenched approximation of QCD, we take the real part of Z(3)-rotated Polyakov loop 
as the Polyakov loop expectation value for the single quark free energy Fq = — Tln(L), because this 
Z(3) branch is chosen when dynamical quarks are present. 
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short distance physics is insensitive to the temperature. Note that, in the case of the 
conventional fixed-A^ approach, because the simulation point is different for each 
T, F 1 (r, T) gets different renormalization depending on T. Therefore, we cannot 
directly compare the results of the bare free energies. In early studies, this in- 
sensitivity was just assumed and used to adjust the constant term of F 1 (r, T) at 
each rPDJUJ More recently, the renormalization factors are estimated by carrying 
out zero-temperature simulations at each finite temperature simulation point in the 
fixed- Nt approach. 35 ' In the present study with the fixed-scale approach, since the 
coupling parameters are common to all T's, we can directly compare bare free ener- 
gies at different T and confirm the expected insensitivity at short distances without 
additional adjustments. 

At large r, we find that i ?1 (r, T) departs from V(r) and eventually becomes flat 
due to the Debye screening. In Figs. [I] and [21 the asymptotic values of F 1 (r, T) at long 
distance are compared with twice the single-quark free energy, 2Fq(T), denoted by 
the arrows on the right side of each plot. We find that i ?1 (r, T) converges to 2Fq(T) 
quite accurately at long distances. This previously was observed also in two-flavor 
QCD with the fixed- Nt approach!^ 1 

Comparing the results of the static-quark free energies in Figs. Q] and [21 we find 
no qualitative differences between quenched and 2 + 1 flavor QCD. In Sec. [H we 
perform a more quantitative comparison in terms of the screening masses. 

Let us now focus on the color-channel dependence by introducing normalized free 
energies V M {r,T) = F M (r,T) - 2F Q (T) (with M = 1,8,6,3*), which are defined 
such that V M = at r — > oo. This is equivalent to define the free energies by 
dividing the right-hand sides of Eqs. (pRjl - tpbrj) by (TrJ?) 2 , because F M (r, T) for all 
channels converges to 2Fq(T) at r — > oo. Results in 2 + 1 flavor QCD are shown in 
Fig. [3] for color-singlet and octet QQ channels (left) and color-antitriplet and sextet 
QQ channels (right). 

We find that the inter-quark interactions, which are Debye-screened at long 
distances, show "attraction" in the color singlet and antitriplet channels, while they 
show "repulsion" in the color octet and sextet channels, at intermediate and short 
distances. These features have been reported also in quenched QCD^ and in two- 
flavor QCD with improved Wilson quarks^ in the fixed-A^ approach. 

§4. Screening masses 

To study screening properties at long distances in each color-channel more quan- 
titatively, we fit the static-quark free energies by the screened Coulomb form 

V M (r,T) = C M ^^-e~ m ^ r , (4-1) 
r 

where a e g (T) and m u (T) are the effective running coupling and the Debye screening 
mass, respectively. Note that V M = at r — > oo with this ansatz, and thus we fit 
the normalized free energies by (|4-ip . 

In the previous studies in quenched!® and two-flavor QCDp^''® the color- 
channel dependence of the static-quark free energies in the high temperature phase 
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Fig. 3. Normalized free energies for color-singlet and octet QQ channels (left) and color-antitriplet 
and sextet QQ channels (right) in 2 + 1 flavor QCD with the improved Wilson quark action. 



turned out to be well absorbed by the kinematical Casimir factors Cm = (X/a=i ^i ' 



C 8 



1 

6' 



C fi 



1 

3' 



2 
3' 



(4-2) 



as inspired from the lowest order perturbation theory. This means that, by fixing 
the coefficient Cm in (|4- 1|) to the value given in (|4-2|) . we can well fit the normalized 
free energies with a e s(T) and m£>(T) which are common to all color channels. With 
our present data, we find that the normalized free energies in both quenched and 
2 + 1 flavor QCD have the same property. Therefore, in the followings, we fix the 
coefficient Cm to the kinematical Casimir factor (j4-2j) to extract a e s(T) and m£>(T). 



4.1. Screening mass in quenched QCD and comparison with fixed-Nt approach 

We first study the Debye screening masses in quenched QCD with the fixed- 
scale approach, and compare them with those obtained in a previous study in the 
fixed- Nt approach. To extract the Debye screening masses we fit the color-singlet 
free energies shown in Fig. Q] with the screened Coulomb form (|4-ip . The fit is 
performed by minimizing x 2 /-^df with the fit range of 1.0 < rT < 1.5, and we obtain 
X 2 /-/Vdf ~ 1-2, except for highest temperature where we obtain xV-^df = 0(100). 
An origin of the large xV-^df at high temperature is that we employ a fixed rage 
in rT for the fits: The number of data within the fit range decreases with increasing 
temperature. Results of mo{T) are shown in Fig. 2] as a function of temperature. 
We also show the results of mi)(T) obtained in the fixed- Nt approach with the tree 
level Symanzik-improved gauge action on iY"i? x N^^ = 32 3 x 4 latticed We take a 
common value of T c for both approaches. We find that the temperature dependence 
of m/j(T) in the fixed-scale approach is quantitatively consistent with that in the 
fixed-iVi approach^ except for the data at the highest temperature T ~ 3.5T C . A 
similar discrepancy between the two approaches at high T has been reported also 
for the equation of stateP^ 

The discrepancy at the highest temperature should be regarded as a consequence 
of the limitation of the fixed-scale approach, since the lattice spacing a may be too 
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Fig. 4. Debye screening mass for the singlet channel in quenched QCD. Open squares are the results 
obtained in the fixed-scale approach. For comparison, results obtained in the fixed- Nt approach 
at Nt = 4^ are also shown by open circles. 



coarse to resolve thermal fluctuations at very high T's in the fixed-scale approach. 
We may estimate the range of applicability of the fixed-scale approach as follows: Let 
us assume that a typical thermal wave length A at finite T is given by A ~ 1 / E where 
E is an average energy of massless particles. We estimate E ~ 3T£(4)/£(3) ~ 2.7T 
for the Bose-Einstein distribution and E ~ 3TC(4)/C(3) x 7/6 ~ 3.15T for the 
Fermi-Dirac distribution. We thus obtain A ~ 1/(3T). On the present quenched 
lattice, A becomes as small as 0.65 fm at T ~ 3.5T C , while our fixed lattice spacing 
in the quenched simulations is a s ~ 0.97 fm. This is probably too coarse to extract 
in-medium properties of thermal matter. 

We should, therefore, take care with the data at temperatures T <; l/(3a). This 
is a potential drawback of the fixed-scale approach around the high temperature 
limit. In the fixed- Nt approach, on the other hand, the lattice becomes finer as 
we increase T, and thus the resolution for thermal fluctuation is kept fixed. In the 
fixed-iVi approach, however, the spatial volume may become quite small at high 
temperatures unless N s is increased simultaneously. 

4.2. Debye screening in 2 + 1 flavor QCD 

We now study the Debye screening in 2 + 1 flavor QCD by performing the 
screened Coulomb fits to the normalized free energies in all color channels. Consult- 
ing the effective screening masses calculated in Appendix [Aj we choose the interval 
0.6 < rT < 1.25 as the fit range. The fit is performed by minimizing x 2 /Nf>f- Re- 
sults of x 2 /Ndf for each color-channel and temperature are summarized in Table ITTTl 
We find that the values of x 2 /Ndf at the highest temperature (T = 697 MeV) are 
quite large in contrast to other temperatures. This is due to a fixed range in rT 
discussed above. In Appendix [Aj we also perform fit with a fixed range in r/a. We 
find that the results of this alternative fit are consistent with those of the present fit, 
and the systematic error estimated by the difference due to choice of the fit range is 
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Fig. 5. Results of effective running coupling a e a (left) and Debye screening mass m_o (right) for 
each color-channel as a function of temperature in 2 + 1 flavor QCD. 



Table III. Results of y 2 /Nuf in 2 + 1 flavor QCD for each color-channel and temperature. The fits 
at T = 232 MeV for M = 8 and 6 are unstable since signals at these parameters are smaller 
than their statistical errors in the fitting interval. 



T [MeV] M = 1 8 6 3* 



232 


1.42 






2.43 


279 


0.80 


1.58 


1.38 


2.12 


348 


1.05 


1.18 


1.82 


2.53 


465 


4.36 


0.85 


1.95 


3.48 


697 


69.36 


8.96 


22.92 


48.77 



at most a few percents even at the highest temperature. 

As the limiting temperature T ~ l/(3a) for the fixed-scale approach mentioned 
in the previous subsection, we have about 930 MeV for our lattice spacing 1/a = 
2.79(1) GeV in 2 + 1 flavor QCD. Our highest temperature 697 MeV is lower than 
that. However, since it is of the same order of this scale and Nt = 4 there is rather 
small, we need to keep in mind possible influence of lattice artifacts around the 
highest temperature. 

Results of a e Q (T) and m p (T) are summarized in Table IIVI and plotted in Fig. [5j 
We see that both a e s(T) and rri£,{T) are well independent on the color channel at 
T > 348 MeV, as implied by the Casimir dominance at high temperatures discussed 
at the beginning of this section. This suggests that spatial correlations in different 
color-channels are dictated by the same effective dynamics in the Coulomb gauge. 

To make a qualitative comparison of the Debye masses between lattice QCD and 
thermal perturbation theory, in the left panel of Fig. [6l we show our Nf = 2 + 1 data 
of moiT) in the color-singlet channel, together with the results of leading order (LO) 
and next-to-leading order (NLO) thermal perturbation theory in massless Nf = 3 
QCD. See Ref. 9) for the LO and NLO formulae we adopt. We find that NLO results 
and the lattice data are qualitatively consistent with each other, though there is 
always a concern about the convergence of thermal perturbation theory. The slight 
deviation at the highest temperature T = 696 MeV may be due to the lattice artifact 
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Table IV. Results of a e ff(T) and moiT) in 2+ 1 flavor QCD. The statistical errors are determined 
by a jackknife method with the bin-size of 20 trajectories. The second parentheses in the color- 
singlet channel show systematic errors due to the fit ranges, as explained in Appendix 1X1 







T [MeV] 


M = 1 


8 


6 


3* 


232 


2.86(181)(37) 






0.22(29) 


279 


1.16(26) (33) 


0.06(2) 


3.25(424) 


0.67(27) 


348 


0.80(2)(0) 


0.87(33) 
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0.63(1)(14) 
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Fig. 6. (Left) The Debye screening mass on the lattice in the color-singlet channel together with 
that calculated in the leading-order (LO) and next-to-leading-order (NLO) perturbation theory 
shown by dashed-black and solid-red lines, respectively. The bottom (top) line expresses a result 
at n = 7rT (3ttT), where fi is the renormalization point. (Right) Flavor dependence of the Debye 
screening masses. We assume the pseudo-critical temperature for 2 + 1-flavor QCD as T c ~ 190 
MeV. 

mentioned above. 

Finally we study the flavor-dependence of the Debye screening masses. In the 
right panel of Fig. [6l we compare the results of m.o(T) in the color-singlet channel 
in 2 + 1 flavor QCD, two-flavor QCD 9 ' and quenched (0-flavor) QCD as a function 
of T/T c , where T c is the (pseudo) critical temperature. The two-flavor results are 
obtained in Ref. ED with improved Wilson quarks in the fixed- Nt approach on a 
16 3 x 4 lattice, where the light quark mass is mpg/my = 0.65. For 2 + 1 flavor QCD, 
we assume T c ~ 190 MeVP^ Note that the two-flavor data are obtained with the 
fixed- Nt approach with the cutoff scales changing from a ~ 0.26 fm (at T ~ T c ) to 
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0.07 fm (at T ~ 4T C ), while the 2 + 1 flavor data in the fixed-scale approach are 
obtained at fixed a ~ 0.07 fm. The two-flavor data may suffer from cutoff effects at 
low temperatures. Consulting the right panel of Fig. [UJ we find that magnitude of 
run in 2 + 1 flavor QCD is about 1.5 times as large as that in quenched QCD, except 
for the data at T ~ 3.5T C in quenched QCD which suffers from lattice artifacts as 
discussed in the previous subsection. (The limiting temperature T ~ l/(3a) for the 
fixed scale approach is about 2.3T C for the quenched data, and about 4.9T C for the 
2 + 1 flavor data.) 

This order of the magnitude for m£>(T) in quenched and 2 + 1 flavor QCD is 
roughly consistent with the prediction of the thermal perturbation theory in which 
the screening mass is predicted to be m^P = a/T+ Nf/6 gT in the leading order 
where g is the running coupling and the quark mass effects are neglected*^!. On the 
other hand, we find that the magnitude of mo in 2 + 1 flavor QCD is almost com- 
parable to that in two-flavor QCD at mpg/my = 0.65,® which is approximately the 
same value with that in 2 + 1 flavor QCD. This suggests that the dynamical strange 
quark does not have large effect for the Debye screening mass in the temperature 
region we investigated. 

§5. Summary 

We have studied static-quark free energies and Debye screening masses in finite 
temperature lattice QCD adopting the fixed-scale approach. Unlike the conventional 
fixed- Nt approach, the fixed-scale approach enables us to unambiguously estimate the 
temperature dependence of physical quantities without varying the spatial volume 
and renormalization factors. On the other hand, the fixed-scale approach may suffer 
from lattice artifacts at T ^ l/(3a) where the thermal wave length of the plasma 
particles become comparable to the lattice spacing. 

Performing simulations in quenched QCD and also in 2 + 1 flavor QCD with 
non-perturbatively 0(a)-improved Wilson quarks in the high temperature phase, we 
found that the static-quark free energies in the color-singlet channel converge to the 
zero temperature static-quark potential at short distances. This confirms the theo- 
retical expectation that short distance physics is insensitive to the temperature. At 
long distances, the free energies approach to twice the single-quark free energies, im- 
plying that the interaction between static quarks is fully screened. We also confirmed 
that the inter-quark forces are "attractive" in the color-singlet QQ and antitriplet 
QQ channels and "repulsive" in the color-octet QQ and sextet QQ channels in the 
case of 2 + 1 flavor QCD. The static-quark free energies in these color-channels are 
well fitted by the screened Coulomb form with the kinetic Casimir factors, which are 
color-channel dependent, combined with color-channel independent effective coupling 
a cs(T) and Debye screening mass m,£,{T). These observations are in accordance with 
previous studies in quenched^ and two-flavor QCD^ in the fixed-iVj approach. Our 

*' A more quantitative comparison may be attempted as follows: Assuming one-loop running 
coupling g~ 2 (/i) = (4-7r)~ 2 (11 — §JV/) In (fi/A) 2 and disregarding a small difference in Ap' we find 
rriD{Nf — 0) : mo(Nf — 2) : mo{Nf = 3) ~ 1 : 1.23 : 1.35 in the massless limit, to be compared 
with the lattice result ~ 1 : 1.5 : 1.5 shown in Fig. [S] (left) at finite quark masses. 
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result of rri£)(T) is qualitatively consistent with the next-to- leading order thermal 
perturbation theory except at the highest temperature corresponding to Nt = 4. 
Values of m£>(T) in 2 + 1 flavor QCD are larger than those in quenched QCD, which 
is again qualitatively consistent with a prediction of the thermal perturbation the- 
ory. We also noted that m£>(T) in 2 + 1 flavor QCD are comparable to those in the 
two-flavor QCD obtained in the fixed- approach at a similar light quark masses,® 
implying that, at the quark masses of the present study, the influence of the dynam- 
ical strange quark is not large in the Debye screening mass in the temperature range 
we investigated. 
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In order to extract the screening mass rrt£), we fit V (r,T) by the screened 
Coulomb form, Eq. (|4-ip . To determine an appropriate fit range, we calculate the 
effective screening mass defined by 



for the color-singlet channel. Figure [7] shows the results of effective screening masses 
at several temperatures as a function of rT. Since we find that the effective masses 
have a plateau at mid range of the distance, we choose the fit range to be 0.6 < 
rT < 1.25. The fit results of the screening masses shown as the lines on each 
figure reproduce the plateau well. Results of xV-^df for each color-channel and 
temperature are summarized in Table IIIII 

Although Fig. [7] suggests that the plateau region in r is mainly determined by 
the thermal length 1/T, we may attempt to fit the screening masses adopting a 
range independent of T. We find that the range 4.5 < r/a < 8.0 also leads to fits 
with acceptable x 2 /-^df- We estimate systematic errors due to a variation of the 
fit range by taking the difference of the fit results between 0.6 < rT < 1.25 and 
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Appendix A 

Effective screening masses and systematic errors 



m<§(T;r) = — log 



1 V x (r) 1 




Ar V 1 (r + Ar) Ar 
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Fig. 7. Effective Debye screening mass at several temperatures as a function of rT. The shaded 
lines present the fit range with the fit results and statistical errors. 



4.5 < r/a < 8.0. The systematic errors thus estimated are given in the second 
parentheses of a e s(T) and m£)(T) in Table IIVI for the color-singlet channel. We 
find that the systematic errors are comparable to or smaller than statistical errors, 
except for the highest temperature. The magnitude the systematic error at the 
highest temperature is, however, at most a few percents. 
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